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THE GENERALIZED MARCUM Q-FUNCTION: AN ORTHOGONAL 

POLYNOMIAL APPROACH 



SZILARD ANDRAS, ARPAD BARICZ, AND YIN SUN 



Abstract. A novel power series representation of the generalized Marcum Q— function of pos- 
itive order involving generalized Laguerre polynomials is presented. The absolute convergence 
of the proposed power series expansion is showed, together with a convergence speed analysis 
by means of truncation error. A brief review of related studies and some numerical results are 
also provided. 

7— I . 

o 

{Sj ■ 1. Introduction 

For v real number let I v be denotes the modified Bessel function [49, p. 77] of the first kind 
r T ^ ■ of order v, defined by 

(1.1) I v {t) ~ ^ it,2)2n+V 



< 



^ n\T(i< + n + 1) ' 

n>0 v ; 



and let b \-t Q u (a, b) be the generalized Marcum Q— function, defined by 



1 f°° t 2 +a 2 

(1.2) Qu(a,b) = — I t u e ~ I v -x{at)dt, 

> . a Jb 

where b > and a, v > 0. Here T stands for the well-known Euler gamma function. When v = 1, 
the function 

*2 . „2 



bh^Q 1 (a,b)= / te 2 I (at)dt 
Jb 

is known in literature as the (first order) Marcum Q— function. The Marcum Q— function and 
its generalization are frequently used in the detection theories for radar systems [27] and wire- 
less communications [12, 13], and have important applications in error performance analysis 
of digital communication problems dealing with partially coherent, differentially coherent, and 
non-coherent detections [38, 40]. Since, the precise computations of the Marcum Q— function 
and generalized Marcum Q— function are quite difficult, in the last few decades several authors 
worked on precise and stable numerical calculation algorithms for the functions. See the papers 
of Dillard [14], Cantrell [7], Cantrell and Ojha [8], Shnidman [34], Helstrom [17], Temme [46] 
and the references therein. Moreover, many tight lower and upper bounds for the Marcum 
Q— function and generalized Marcum Q— function were proposed as simpler alternative evalu- 
ating methods or intermediate results for further integrations. See, for example, the papers of 
Simon [35], Chiani [10], Simon and Alouini [37], Annamalai and Tellambura [1], Corazza and 
Ferrari [11], Li and Kam [22], Baricz [4], Baricz and Sun [5, 6], Kapinas et al. [19], Sun et al. 
[41], Li et al. [23] and the references therein. In this field, the order v is usually the number of 
independent samples of the output of a square-law detector, and hence in most of the papers 
the authors deduce lower and upper bounds for the generalized Marcum Q— function with or- 
der v integer. On the other hand, based on the papers [8, 27, 34] there are introduced in the 
Matlab 6.5 software the Marcum Q— function and positive integer order generalized Marcum 
Q— function 1 : marcumq(a,b) computes the value of the first order Marcum Q— function Q±(a, b) 
and marcumq(a,b ,m) computes the value of the mth order generalized Marcum Q— function 
Q m (a, 6), defined by (1.2), where m is a positive integer. However, in some important appli- 
cations, the order v > of the generalized Marcum Q— function is not necessarily an integer 
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^See http://www.mathworks.com/access/helpdesk/help/toolbox/signal/marcumq.html for more details. 
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number: The generalized Marcum Q— function is the complementary cumulative distribution 
function or reliability function of the non-central chi distribution with 2u degrees of freedom 
[18, 39, 41]. Moreover, real order generalized Marcum Q— function has been used to charac- 
terize small-scale channel fading distributions with line-of-sight channel components [24, 50] or 
cross-channel correlations [2, 3, 19, 20, 38, 44, 45]. 

In this paper, we present a novel generalized Laguerre polynomial series representation of the 
generalized Marcum Q— function, which extends the result of the first order Marcum Q— function 
in Pent's paper [32] to the case of the generalized Marcum Q— function with real order v > 0. We 
further show the absolute convergence of the proposed power series expansion, together with a 
convergence speed analysis by means of truncation error. A brief review of related studies in the 
literature is provided, which may assist the readers to get a more complete vision of this area. 
Finally, some numerical results are provided as a complementary of these theoretical analysis. 

2. The generalized Marcum Q— function via Laguerre polynomials 

2.1. Novel series representation of the generalized Marcum Q— function. We start with 
the following well-known formula [43, p. 102] 

(2.1) £ j^ Z ~t = r(a + l)e 2 (x,)-f J a (2Vx~z), 

n>0 L n (0) 

where x,z£l and a > — 1. Here J a stands for the Bessel function of the first kind of order a, 
is the generalized Laguerre polynomial of degree n and order a, defined explicitly as 

T^)(r\ = ( p -^ T n+a)(n) = j> + a + 1) {- X f 

n { } n! [ > ^T(k + a + l)T(n-k + l) k\ ' 

k=0 

Changing in (2.1) z with — z and taking into account the formula I u {x) = \~ u J u (ix) we obtain 
that [26] 

L ( n\x) (-l) n Z n 
i>0 \}J) 



(2.2) £ = r(a + l)e-*M-f 7.(2^). 



Now, if we use 



lW(o)- r (n + a + i) 



r(a + l)r(n + l)' 

and replace x with a and a with v — 1, respectively, (2.2) can be rewritten as 
(2.3) (£) V e- z -%_!(2^) = e' a ^(-l) n L J: ^H z**"- 1 



T(v + n) 

n>0 K ' 



which holds for all a, v > and z > 0. 

Now, consider the following formula [46, 47] 

T 1 1 i Tl, POO v — 1 

(2.4) Q u (V2-a, V2b) = e - £ ^ + /T = / (*)~ ir e— 7,_ 1 (2 v ^)d,, 

n>0 v ' Jb 

where a, v > and b > 0. We note that the function 6 h-> Q v (\fa., Vb), defined by 

Q v (y/a,Vb) = ^ y 2 e-^/^itVo^dz, 

is in fact the survival function (or the complementary of the cumulative distribution function 
with respect to unity) of the non-central chi-square distribution with 2v degrees of freedom and 
non-centrality parameter a. With other words, for all a, v > and b > we have 

(2.5) Qu(V^,Vb) = l--J^ (£) 2 e-^i^V^ds. 
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See [39] for more details. Combining (2.3) with (2.5) we obtain 



Q„(V2^,V2b) = l- / (-) 2 e-^I^^^dz 
Jo 

rb T {y ~ V> (n\ 

= 1-1 e -y ( _i)«^ \?l z n+»-i dz 

Jo f^n T{u + n) 



4o r ( , + n+i) 

where in (a) the integration and summation can be interchanged, because the series on the 
right-hand side of (2.3) is uniformly convergent for < z < b. For more details see the last 
paragraph of Section 2.2. After some simple manipulation, we obtain a new formula of the 
generalized Marcum Q— function, i.e., 

Llr 1 ' (€) /«\«+- 



(2.6) W) = 1 _ E( _ ire -#_ilZ(y 

n>0 

valid for all a, v > and 6 > 0. 

In order to simplify the numerical evaluation of the series (2.6), we consider the expression 

m -" nL ^ ){a) 



I> + n+l)' 

which satisfies the recurrence relation 

d / j,\ (2n + z/-a)6 

P„ )n+ i(a,&) =7 — —— — ■ —P^ n {a,b) 

(n + l)[v + n + 1) 

(n + i/ - 1)6 2 
" ( n + + „)(„ + n + i) P »,n-i(a, b) 

for all a, v > 0, & > and n € {1, 2, 3, . . . }, with the initial conditions 

^°M ) = i>TI) ^d P v>l(a> &) = i^g. 

Here, the recurrence relation for P U:n (a, b) were obtained from the recurrence relation [43, p. 
101] 

(n + l)4^(x) = (2n + a + 1 - x)4 a )(s) - (n + <*)4-i 0*0 
and the initial conditions from 

4 Q) (x) = l and 4 Q) (x) = -x + q + 1. 
With the help of the expression P u , n (a, b), (2.6) can be easily rewritten as 

(2-7) Q„(a, &) = 1 - J>"* (|) " P W>B - |) . 

n>0 ^ ' ^ ' 
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2.2. Convergence analysis of the new series representation. We note that for a > 0, 
v > 1 and 6 > the absolute convergence of the series in (2.6) or (2.7) can be shown easily by 
using the following inequalities 



or 



£(-ire-i 



n »-V(a*\ ... 



n>0 



2\ n +v 



a x — ^ 



r> + n + 1) V 2 

1 / fo 2X«+. 



n>0 



r(z/ + n + 1) V 2 



a x — ^ 

<e-r£ 



1 



n>0 



2\ 



< e * 



e 4 



i» 



i» 



r(i/ + n+l) V 2 

1,2 \ "-I 

T E 



1 



L („-l) 



r(i/ + n) «i 

— -e 4 

n!r(i/) 

2 \ "+ 1 



n>0 

e 2 



(n + 1)! V 2 



2 L^r^i 



n>0 



r(i/ + n + 1) V 2 



- 6 2 ^ I> + n + 1) V 2 

n>0 



2 \ n+f 



< e 2 



E 

n>0 



1 



T{v + n+ 1) V 2 



L (-i) 
n 

2\ n+f T-l/ 

1 (z/ + n 



2b 



2\ 



< e -v " V 



e 4 



or 

262 
a 2 



1 



^ (n + 1)! 

n>0 V ; 



n! 

2 



e 2 



e 4 — 



1-1/ 



which contain the known inequalities of Szego [43] for generalized Laguerre polynomials 



\L a (x)\ < r(a + re + 1) c f 
|L " (X)| - n!r(a + l) 6 



and of Love [25] 



\K(x)\ < 



T(a + n + 1) fx\ - a x 



where in both of the inequalities a > 0, x > and n £ {0, 1,2,... }. 
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Moreover, for a > 0, < v < 1 and b > the absolute convergence of the series in (2.6) or 
(2.7) can be shown in a similar manner by using the following inequality 



E(- 1 )"^ 

n>0 



r("-l) 
„2 Li n 



< e 



T(u + n+ 1 
1 



n+v 



n>0 



n+v 



< e 



E 

n>0 

£ 

n>0 

£ 



1 



r(i/ + n + l) 



1 



2J 



L ("-i) 

r(z/ + n 



1 



z^ + n VT(^ + n) n!r(i/) 



n!r(i/ 

2 \ n+v 



e 4 



2e" 



n>0 



n+v 



2 ] 62 ' 



which contains the classical inequality of Szego [43] for generalized Laguerre polynomials 

T(a + n + 1) 



|L«(x)| < 2 



n!r(a + l) 



e2 



where — 1 < a < 0, x > and n£ {0,1,2,... }. In addition here we used the fact that for all 
fixed n G {1, 2, 3, ... } the function 



1 



1 



i/ + n \T{u + n) n\T(v) 



which maps into 2/n!, is decreasing on (0,1] and consequently for all n G {0,1,2,. 
< v < 1 we have 



} and 



1 



1 



v + n \T(u + n) n\T(y) 



< 



We note that other sharper uniform bounds for generalized Laguerre polynomials can be found 
in the papers of Love [25], Lewandowski and Szynal [21], Michalska and Szynal [28], Pogany and 
Srivastava [33]. See also the references therein. 

Finally, note that by using the above uniform bounds for the generalized Laguerre polynomials 
the uniform convergence of the series on the right-hand side of (2.3) can be shown easily for 
< z < b. This is important because in order to obtain (2.6) we have used tacitly that the 
series on the right-hand side of (2.3) is uniformly convergent and then we can interchange the 
integration with summation. For example, if we use the above Szego's uniform bound, then for 
all n G {0, 1, 2, . . . }, a > 0, v > 1 and < z < b we have 



V ; r(i/ + n) 



< 



T(u) n! 



By the ratio test the series e b = y^ n>0 &™ /n\ is convergent and thus in view of the Weierstrass 
M-test the original series on the right-hand side of (2.3) converges uniformly for all < z < b. 



2.3. Truncation error analysis. For practical evaluations of our power series expansion, we 
need to approximate the generalized Marcum Q— function Q u (a, b) by the first uq G {1, 2, 3, ... } 
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terms of (2.6), i.e., 



n 



Q v {a,b) = l-J2(-V n e- S * 



n=0 



r("-l) (a 2 

I> + n + 1) 



52 \ 



We note that the absolute value of the truncation error 

e t = Q u (a,b)-Q u (a,b)= £ (-lf+V^ 



r("-l) 
2 -l^n 



n>no+l 



I> + n + l) V 2 



2 \ n+u 



can be upper bounded by using the upper bounds for the generalized Laguerre polynomials as 
in subsection 2.2. More precisely, by using the same argument as in subsection 2.2 and SewelPs 
inequality [29, p. 266] 



Ex xe x 
H<—> n G {l,2,3,...}, x>0, 



k=0 



xe 
n 



we can deduce the following: if a > 0, b > and v > 1, then 



M < e 4 



r» 



2\ v —l 



or 



(2b 

N < e 4 — 



2 

2 \ 



6_ 

e ^ 



, "0+1 -1 



n ! ( 9 



n=0 



< 



e 2 



1 



n + 1 I» 



ft2\ 



e 2 



no+1 1 



-E-(- " 

n=0 



< 



e 2 



,2 2 

-"^6 2 / 26^ 1,-1 



n + l 2 



Similarly, it can be shown that if a > 0, b > and < v < 1, then the absolute value of the 
truncation error is upper bounded as follows 



2\ " 



n 



e 2 



-E 



1 



n=0 



n! 



2\ n 



< 



2e~ 



n 



2\ "+1 



Observe that the above upper bounds of the absolute value of the truncation error converge to 
zero at a speed of 1/tiq. In practice, we can use these upper bounds to decide the number of 
terms, i.e. no, for achieving a pre-determined accuracy. 

2.4. A brief review of related studies. As far as we know the formula (2.6), or its equivalent 
form (2.7), is new. However, if we choose v = 1 in (2.7), then we reobtain the main result of 
Pent [32] 



Qi(a,b) 



Ea / U U 



n>0 



where 



P n {a,b) = P hn (a,b) 



2 ' 2 
b n L n (a) 



(n + 1)!' 

which for all a > 0, b > and nG {1,2,3,...} satisfies the recurrence relation 



P n+1 (a,b) 



with the initial conditions 



Here L n = Ln^ is the classical Laguerre polynomial of degree n. 



(2n + l -a)b 
(n + l)(n + 2) 



P n (a,b) 



n 



b 2 



P (a,b) = l and Pi (a, b) = 



(n + l) 2 (n + 2) 
(1 - a)b 



P n -i(a,b) 
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It should be mentioned here that another type of Laguerre expansions for the Marcum 
Q— function was proposed in 1977 by Gideon and Gurland [16], which involves the lower in- 
complete gamma function. This type of Laguerre expansions requires to use a complementary 
result of (2.2), i.e. 

(2.8) £ L l) { * n ~ l)n * n = r (« + l)e- x (xz)-fl a (2VTz). 

n>0 {V) 

Now by some simple manipulation we obtain 

(2.9) (£) 2 r^. 1 (2^ = rV«^±i,t 1) W- 

a n>0 { + > 

which is equivalent to Tiku's result [48], available also as equation (29.11) in the book [18]. By 
integrating (2.9) in z and by using the differentiation formula [26] 

Tz r 

where n £ {1,2,3, ...}, a > — 1 and z 6 R, we can obtain another generalized Laguerre polyno- 
mial series expansion of the generalized Marcum Q— function 



z a+l e- z L { ^\z)] = nz a e- z L^(z), 



o", 



which in turn implies that 



-r(u) T {^2) 2J !) e [ 2 ) nr( , + n)U 



where 7(-, •) is the lower incomplete gamma function, defined by 

7 (a,x) = I t^e^dt. 
Jo 

Here we used that 

(2.11) T(a,x) = T(a)- 1 (a,x), 

and 

SS°" (a ' 6) = i^) r ("'7 

Some other Laguerre expansions for the Marcum Q— function are provided in Gideon and Gur- 
land's paper [16], available also as equation (29.13) of [18]. Moreover, a new unified Laguerre 
polynomial-series-based distribution of small-scale fading envelope and power was proposed 
recently by Chai and Tjhung [9], which covers a wide range of small-scale fading distribu- 
tions in wireless communications. Many known Laguerre polynomial-series-based probabil- 
ity density functions and cumulative distribution functions of small-scale fading distributions 
are provided, which include the multiple-waves-plus-diffuse-power fading, non-central chi and 
chi-square, Nakagami-m, Rician (Nakagami-n) , Nakagami-g (Hoyt), Rayleigh, Weibull, Stacy, 
gamma, Erlang and exponential distributions as special cases. See also [42], which contains 
some corrections of formulas deduced in [9]. In particular, (2.10) is a special case of the unified 
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cumulative distribution function given in corrected form in [42]. We note that the expression 
of (2.10) and the unified cumulative distribution function in [42] are quite different from our 
main result (2.6) or (2.7). This is because they are based on two different Laguerre polynomial 
expansions of the modified Bessel function of the first kind I u given in (2.2) and (2.8). Therefore, 
these Laguerre polynomial expansions are expanded over different variables of the generalized 
Marcum Q— function. Finally, we note that since Nakagami's work [30] the Laguerre polynomial 
series expansions of various probability density functions have been derived. We refer to the 
papers of Esposito and Wilson [15], Yu et al. [51], Chai and Tjhung [9] and to the references 
therein. 

Finally, by using the infinite series representation of the modified Bessel function of the first 
kind (1.1) and the formula 



f 

J a 



n-i ( m + 1 a 2 

2 1 



2 ' 2 

where T(-, •) is the upper incomplete gamma function, defined by 

poo 

F(a,x) = / f-V'dt, 

J X 

we easily obtain that 

tW x ^ (at) 2n+u - ] 



E n 2n poo 2 
^ / „-V/2n+i/-l H/ 
>0 2^+^n\T(u + n) J b 



si 

2 



n>0 



T{u + n) 



(2.12) =!-E 



2 / a 2\ n 7 [v + n 



e 2 

n>0 



r(i/ + n) 



We note that (2.12) is usually called the canonical representation of the vih order generalized 
Marcum Q— function. Recently, Annamalai and Tellambura [2] (see also [3]) claimed that the 
series representation (2.12) is new, however it appears already in 1993 in the paper of Temme 
[46]. See also Temme's book [47] and Patnaik's [31] result from 1949, which can be found also 
as equation (29.2) in the book [18]. Interestingly, our novel series representation (2.7) for the 
generalized Marcum Q— function resembles to the series representation (2.12). 

2.5. Numerical results. We now consider some numerical aspects of our generalized Laguerre 
polynomial expansions (2.6) or (2.7). In practice, we usually need to compute the detection 
probability for different values of b with fixed a to decide a proper detection threshold. Since 
the generalized Laguerre polynomial in (2.6) is determined by only a, we can save computation 
time by storing the values of the generalized Laguerre polynomials for computing the generalized 
Marcum Q— function with different values of b. 

The following tables contain some values of the generalized Marcum Q— function calculated 
using (2.7) and using the Matlab marcumq function. For the considered choices of a and b, 
the numerical value of (2.7) is exactly the same with that of the Matlab marcumq function, 
if v = 1,3,5 is integer. When u = 7.7, the Matlab marcumq function does not work, and 
the numerical value of (2.7) is provided in the tables. Finally, we note that more accurate 
intermediate terms are required for larger a and b. 



a = 0.2, b = 0.6 


v = 1 


v — 3 


v = 5 v = 7.7 


(2.7) 


0.838249985438908 


0.999166310455636 


0.999998670306184 0.999999999927717 


marcumq 


0.838249985438908 


0.999166310455636 


0.999998670306184 
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a = 1.2,6= 1.6 


v = 1 


v = 3 


v = 5 


v = 7.7 


(2.7) 


0.501536568390858 


0.916936068900377 


0.994346394491553 


0.999944937223540 


marcumq 


0.501536568390858 


0.916936068900377 


0.994346394491553 






a = 2.2,6 = 2.6 


v = \ 


v = 3 


v = 5 


v = 7.7 


(2.7) 


0.426794627821735 


0.746459898209090 


0.929671935077756 


0.993735633182201 


marcumq 


0.426794627821735 


0.746459898209090 


0.929671935077756 
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